AlllQM  053343 


NIST 

PUBLICATIONS 


Fracture  Mechanics 
Evaluation  of  Railroad 
Tank  Cars  Containing 
Circumferential  Cracks 


— QC 

100 

.056 

#5179 

1995 


NISTIR  5179 


Report  No.  27 


George  E.  Hicho 
Akram  Zahoor 
Richard  J.  Fields 
Roland  deWit 

U.S.  DEPARTMENT  OF  COMMERCE 
Technology  Administration 
National  Institute  of  Standards 
and  Technology 
Gaithersburg.  MD  20899 


NIST 


Report  No.  27 


Fracture  Mechanics 
Evaluation  of  Railroad 
Tank  Cars  Containing 
Circumferential  Cracks 


George  E.  Hicho 
Akram  Zahoor 
Richard  J.  Fields 
Roland  deWit 

U.S.  DEPARTMENT  OF  COMMERCE 
Technology  Administration 
National  Institute  of  Standards 
and  Technology 
Gaithersburg,  MD  20899 


April  1993 


U.S.  DEPARTMENT  OF  COMMERCE 
Ronald  H.  Brown,  Secretary 

NATIONAL  INSTITUTE  OF  STANDARDS 
AND  TECHNOLOGY 
Arati  Prabhakar,  Director 


^ ^ .£#W^ 


fat'b'.U  .’^. 

, 'MV. 

' 'fv'  :'■'■■  ^ U 

v 'U  ’*■'  "inv'ip-i*  tsficf*''W 

.-^';ti'vjf  ifi;>='''' 


t \''0i‘,  !k'(5A 


■■'5  .^iiiii'!Ai^‘vil^  "^Jfi'nyi 


Report  No.  27 


FRACTURE  MECHANICS  EVALUATION  OF  RAILROAD  TANK  CARS 
CONTAINING  POSTULATED  CIRCUMFERENTIAL  CRACKS 

George  E.  Hicho,  Akram  Zahoor^ , 

Richard  J.  Fields,  and  Roland  deWit 


Mechanical  Properties  and  Performance  Group 
Materials  Science  and  Engineering  Laboratory 
National  Institute  of  Standards  and  Technology 
Metallurgy  Division 
Gaithersburg,  MD  20899 


ABSTRACT 


Fracture  mechanics  analyses  using  the  J - integral  as  fracture  parameter  were 
performed  on  railroad  tank  cars  made  from  normalized  AAR  TC128  grade  B steel. 
Circumferential  throughwall  cracks  in  the  tank  car  shell  region  were  postulated 
to  determine  the  critical  crack  size  for  axial  tension  loadings  anticipated  in 
service.  Five  loading  cases  were  considered.  These  were:  (1)  coupler  impact, 
(2)  test  pressure,  (3)  start- to-discharge  pressure,  (4)  bursting  pressure,  and 
(5)  vapor  tight  minimum  pressure.  The  analyses  were  performed  for  two  service 
temperatures,  -40  °C  (-40  °F)  and  22  °C  (72  °F) . The  bounding  critical  crack 
length  is  determined  to  be  approximately  30  times  the  tank  car  shell  wall  thick- 
ness for  the  start- to-discharge  event.  The  critical  crack  lengths  for  the 
coupler  impact  force  are  approximately  1.7  times  larger  than  those  for  the 
discharge  pressure.  The  critical  crack  length  is  determined  to  be  the  largest 
for  the  impact  force  case.  The  results  also  show  that  even  at  burst  pressure, 
small  throughwall  cracks  (four  times  tank  shell  thickness)  can  be  tolerated. 

INTRODUCTION 


Since  January  1,  1989,  the  normalized  AAR  TC128  grade  B steel  is  to  be  used  for 
railroad  tank  cars  that  may  carry  hazardous  commodities.  Prior  to  this  date,  the 
TC128  grade  B steel  used  in  tank  cars  was  normally  supplied  in  the  as  rolled  con- 
dition. Previous  research  (1,  2,  3,  and  4)  showed  that  the  mechanical  proper- 
ties, in  particular  the  impact  strength  at  low  test  temperatures  for  the  steel 
in  the  as  rolled  condition  was  markedly  lower  than  normalized  TC128  grade  B 
steel.  Fracture  toughness  tests  conducted  on  these  steels  showed  that  the 
normalized  steel  resists  crack  initiation  better  than  the  as  rolled  steel.  As 
part  of  our  present  program  with  the  Federal  Railroad  Administration,  research 
was  initiated  using  data  generated  at  NIST,  to  assess  the  structural  integrity 
of  the  tank  car  shell  containing  postulated  circumferential  throughwall  cracks 
under  axial  tension  loading.  The  latest  fracture  mechanics  analysis  methods  were 
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used  to  determine  the  critical  crack  size  for  service  conditions  at  -40  °C.- 
(-40  "F)  and  22  “C  (72  °F)  . The  fracture  mechanics  methods  and  analyses 
previously  used  for  the  tank  car  evaluation  were  based  on  linear  elastic  fracture 
mechanics  (LEFM) . Further,  the  solutions  used  were  for  the  simplified  case  of 
a flaw  in  a flat  plate.  Consequently,  the  effects  of  tank  shell  geometry  and 
flaw  orientation  were  ignored  in  the  analysis.  The  analyses  performed  in  this 
report  use  appropriate  J- integral  elastic -plastic  fracture  mechanics  solutions 
for  the  tank  car  geometry,  thus  allowing  accurate  evaluation.  The  J based 
analyses  are  superior  to  those  based  on  LEFM  concepts  because  it  accounts  for 
crack  tip  plasticity  consistent  with  the  material  stress  - strain  curve  and  allows 
assessment  of  the  stability  of  crack  extension . The  report  also  discusses  the 
American  Association  of  Railroads  (AAR)  analysis  method  previously  used  for 
assessing  tank  car  structural  integrity  and  compares  with  the  latest  technology 
used  in  this  report.  Areas  for  improvements  are  discussed. 

ELASTIC -PLASTIC  FRACTURE  MECHANICS  METHOD 


The  fracture  analysis  method  used  in  this  report  is  based  on  the  elastic -plastic 
J integral  and  associated  tearing  modulus  crack  stability  criterion.  Ductile 
fracture  analysis  consists  of  two  components:  crack  driving  potential  (J)  and  the 
material's  resistance  to  crack  extension.  The  crack  driving  potential  is  com- 
puted from  a formula  that  depends  on  the  applied  load,  component  geometry, 
service  temperature,  and  stress  - strain  properties  of  the  material.  The  crack 
extension  resistance  is  measured  from  laboratory  specimens.  This  latter  com- 
ponent is  characterized  by  the  J-resistance  curve  of  the  material.  The  crack 
driving  potential  is  then  compared  to  the  material  resistance  to  assess  the 
component's  fitness  for  service. 

The  Crack  Driving  Potential  (J): 

The  fracture  mechanics  evaluation  of  the  railroad  tank  car  shell  was  performed 
using  the  J integral  elastic -plastic  fracture  parameter.  In  general,  the  J- 
integral  solution  is  expressed  as  (5): 

J = Je  + Jp  (1) 


The  first  term  on  the  right  hand  side  of  Eq . (1)  is  the  elastic  J modified  for 
small  scale  yielding.  The  elastic  J can  be  computed  from  the  mode  I stress 
intensity  factor,  Kj , as  Kj^/E  where  E is  the  elastic  modulus  of  the  material. 
The  second  term  is  the  plastic  component  of  the  J- integral,  which  depends  upon 
crack  size,  component  geometry,  and  plastic  characteristics  of  the  stress  strain 
curve . 

Limited  elastic  and  elastic -plastic  J solutions  for  throughwall  cracks  in 
cylindrical  geometries  were  developed  originally  in  Refs.  6,  7.  These  solutions 
were  later  expanded  to  covei’  a wide  range  of  crack  sizes  and  are  given  in  Refs. 
5,  8.  In  general,  solutions  for  shell  geometries  with  a large  radius- to- 
thickness  ratio  (R/t  greater  than  20)  are  available  only  for  linear  elastic 
fracture  analysis,  and  these  are  limited  to  a few  crack  lengths.  Refs.  9-11. 
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Material  Fracture  Resistance.  The  J -Resistance  Curve: 


The  J-resistance  curve  is  most  commonly  generated  from  compact  tension,  C(T), 
specimens  following  the  procedures  set  forth  in  ASTM  El  152-87  (12).  The  speci- 
men size  is  selected  so  that  it  either  has  the  same  thickness  as  the  component 
thickness  or  has  a thickness  that  would  produce  a conservative  estimate  of  the 
material's  resistance  to  crack  extension.  The  direction  of  crack  extension  in 
the  specimen  is  selected  to  match  that  of  the  flaw  in  the  application. 

While  J-R  curves  can  be  generated  for  the  appropriate  direction  of  crack 
extension  in  the  application,  the  data  are  usually  limited  to  small  crack 
extension  due  to  the  specimen  size  limitation.  A standard  ASTM  compact  tension, 
IT  C(T),  specimen  produces  only  .25  to  .51  cm  ( 0.1  to  0.2  inch)  of  valid  crack 
extension.  When  it  is  not  possible  to  fabricate  a standard  thickness  specimen 
due  to  smaller  component  thickness,  it  is  preferred  that  the  specimen  has 
approximately  the  same  thickness  as  in  the  application. 

Stability  Limit.  The  J-Integral  Tearing  Modulus: 

The  analyses  performed  in  this  report  are  based  on  the  J- integral  tearing  modulus 
(J-T)  method.  The  analysis  method  compares  the  crack  driving  potential  (J)  to 
the  material  resistance  (J-R  curve).  During  stable  crack  growth  these  two 
parameters  are  equal  until  a stability  limit  is  reached.  It  may  be  shown  that 
this  limit  is  given  by  tearing  modulus  related  to  the  slope  of  the  J-resistance 
curve . 

The  J-T  analysis  used  here  was  developed  in  1977  (13).  Initially,  it  was  used 
to  demonstrate  the  concept  of  ductile  tearing  instability  to  specimen  geometries, 
which  included  single  edge  notched  bend,  C(T),  and  center-cracked  tension  speci- 
mens (14)  . This  analysis  was  later  applied  to  cracked  pipes  and  the  analysis  was 
based  on  perfectly  plastic  assumption.  Subsequent  work  (15-19)  removed  this 
assumption  and  developed  the  J and  T solutions  for  throughwall  cracks,  part 
throughwall  cracks,  and  compound  cracks  under  bending,  axial  tension,  or 
torsional  loading. 

The  J-T  procedure  for  determining  the  load  at  incipient  crack  instability  is 
described  below.  Suppose  that  the  maximum  bending  moment  needs  to  be  computed 
for  a specified  circumferential  throughwall  crack  length  in  a cylinder.  The 
first  step  in  performing  this  analysis  is  to  determine  the  applied  J- integral  as 
a function  of  the  applied  loads  for  the  specified  cylinder  geometry  and  crack 
size.  The  applied  J is  calculated  as  a function  of  bending  moment  using  Ref.  5 
solution  as  shown  schematically  in  Figure  la.  Next,  applied  J also  is  computed 
as  a function  of  tearing  modulus  (using  Ref.  5 solutions)  for  the  specified  crack 
length  and  a range  of  bending  moments.  The  applied  J and  T pairs  computed  at 
each  bending  moment  are  then  plotted  to  establish  the  applied  J-T  line  shown 
schematically  in  Figure  lb. 

A J-T  curve  representing  the  material  crack  extension  also  is  determined  and 
plotted  on  Figure  lb.  The  material  J-T  curve  is  developed  from  an  experimental 
J versus  crack  extension  (J-R)  curve,  which  can  be  obtained  from  compact  tension 
specimens  or  specimens  with  same  geometry  as  the  component.  The  value  of 
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material  T at  any  specified  value  of  J is  determined  from  the  slope  of  the.J-R 
curve  as ; 


T = (Aj/Aa)  • (E/a^2)  (2) 

where  E and  are  the  elastic  modulus  and  material  flow  stress,  respectively. 
The  flow  stress  is  a reference  quantity  usually  defined  as  one -half  the  sum  of 
yield  and  ultimate  tensile  strengths.  The  AJ/Aa  is  the  slope  of  the  J-R  curve 
at  any  specified  value  of  J. 

The  intersection  of  the  applied  and  material  J-T  curves  illustrated  in  Figure 
1(b)  defines  the  J and  T values  corresponding  to  crack  instability.  The  maximum 
bending  moment  is  determined  by  entering  Figure  la  with  the  J at  crack 
instability  from  Figure  lb  and  determining  the  corresponding  value  of  bending 
moment.  This  procedure  is  repeated  using  current  crack  length  (original  crack 
length  plus  calculated  crack  extension) , and  a revised  bending  moment  is 
obtained. 

A similar  procedure  can  be  used  to  determine  the  allowable  crack  length  for  any 
specified  load  when  Figure  la  is  replaced  by  a plot  of  J versus  crack  length  for 
the  specified  load.  The  procedure  is  the  same  for  part- throughwall  cracks,  any 
other  flaw  geometry,  or  loading  type.  Further  details  can  be  found  in  Ductile 
Fracture  Handbook,  Ref.  5,  which  also  contains  a sample  problem  and  the  step-by- 
step  solution. 

ASSESSMENT  OF  FRACTURE  ANALYSES  METHODS 


AAR  Fracture  Analyses  Method : 

This  section  describes  the  fracture  analysis  method  developed  by  Fellini  for  the 
Association  of  American  Railroads  (AAR)  (20,21).  The  method  is  based  on  LEFM 
principles  and  defines  procedures  for  both  the  crack  initiation  and  crack  arrest 
analysis.  The  method  and  procedures  have  been  defined  so  that,  with  certain 
assumptions,  it  can  be  used  as  a design  reference  graph.  It  was  developed  for 
application  to  low  and  intermediate  steels  having  less  than  482  MPa  (70  ksi) 
yield  strength. 

The  AAR  method  requires  the  following  basic  information  to  determine  the  critical 
load  or  crack  size. 

- Kj  curve  for  the  dynamic  loading  of  interest,  e.g. , curve  for  the  material 
NDT  temperature  or  NDT  temperature  band  for  the  material 
Kjc  versus  temperature  for  the  material 
Material  tensile  properties 

Stress  intensity  factor  .solution  for  the  flaw  geometry 
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Crack  Arrest  Diagram: 


The  procedure  for  developing  the  crack  arrest  diagram  is  described  below.  The 
procedure  requires  definition  of  two  points:  L and  YC.  These  are  defined  on  the 
Kj^  curve.  Figure  2 (21)  shows  the  Kj^  curve,  which  is  a plot  of  Kj^  versus 
temperature  relative  to  NDT . The  abscissa  of  this  figure  can  be  expressed 
alternately  as  temperature  with  NDT  point  marked  to  distinguish  brittle  to 
ductile  transition.  The  relative  temperature  scale  is  preferred  for  graphical 
methods.  The  NDT  temperature  is  established  by  the  drop-weight  test.  The  test 
is  a fixed  crack  size  test  and  is  conducted  by  dynamically  loading  the  specimen 
to  yield  stress  level.  The  Kj^  curve  has  been  widely  used  in  industry,  and  ASME 
Code  procedures  for  reactor  vessels  are  based  on  such  a curve.  The  point  L is 
defined  by  the  plane  strain  limit  for  specified  section  thickness,  and  represents 
the  largest  value  of  Kj^  allowed  in  a valid  test.  The  point  L takes  on  higher 
values  for  higher  section  thickness.  It  defines  the  point  beyond  which  LEFM 
methods  do  not  apply  because  elastic -plastic  fracture  regime  begins  from  this 
point  to  higher  Kj^  values. 

The  Kj^  value  at  point  L is  calculated  from  the  following  formula: 

^id  = (3) 


where  B is  the  section  thickness,  a , is  estimated  to  be  o,,  -i-  158  in  MPa  units 

yd  y 

(23  in  ksi  units),  where  is  0.2  percent  offset  yield  strength  (20). 

The  YC  point  refers  to  the  yield  criterion  value  of  Kj^  for  specified  section 
thickness.  It  signifies  a condition  of  crack  arrest  for  applied  stress  equal  to 
yield  stress.  The  YC  point  is  relevant  to  structures  subjected  to  accidental 
overloading.  This  point  is  clearly  well  into  the  elastic -plastic  fracture 
region,  and  an  elastic -plastic  fracture  parameter  should  be  used  to  define  the 
YC  point.  In  the  AAR  method  (20),  this  point  is  defined  by  assuming  LEFM 
behavior : 


K 


Id 


(A) 


where  B and  have  the  same  definition  as  defined  for  Eq . (3).  An  examination 
of  Eqs.  (3)  and  (4)  indicates  that  the  YC  and  L values  are  related  to  each  other 
by  a constant  factor  (72 . 5 ) regardless  of  the  section  thickness,  NDT  temperature, 
and  material  properties. 

Figure  3 (20)  shows  the  L-YC  lines  on  a plot  of  applied  stress  to  yield  stress 
ratio  (cr/Oy)  versus  temperature  relative  to  NDT.  The  Kj^  value  at  L and  YC 
points  is  calculated  by  using  Eqs.  (3)  and  (4)  for  specified  section  thickness. 
These  values  are  then  used  with  the  Kj^  curve  of  Figure  2 , giving  the  reference 
temperatures  for  the  L and  YC  points.  These  points  are  shown  on  Figure  3,  where 
the  value  a/a^  for  the  L point  is  set  at  0.2  based  on  the  experimental  data. 
This  is  an  approximation  since  there  is  some  variation  in  this  value  for  dif- 
ferent materials  and  section  thickness.  The  L-YC  curve  is  assumed  as  straight 
line,  although  a curve  can  be  developed  using  LEFM  or  elastic -plastic  fracture 
mechanics  (EPFM)  methods.  The  L-YC  lines  of  Figure  3 form  the  basis  for  the 
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crack  arrest  diagram.  The  crack  arrest  line  for  a 2.54  cm  (1-inch)  thick  section 
is  shown  in  Figure  4 (21)  . It  is  assumed  that  this  curve  applies  also  to  section 
thickness  less  than  2.54  cm,  but  the  justification  is  lacking. 

Crack  Initiation  Diagram : 

The  crack  initiation  curves  are  developed  using  fracture  toughness  at  the 
initiation  of  crack  extension,  K|^  as  the  basis.  It  is  more  appropriate  to  use 
Kj.  instead  of  Kj^  for  thinner  sections.  The  Kj  calculations  are  performed  for 
a specified  flaw  size  with  applied  load  (or  stress)  as  a parameter.  The  calcu- 
lated Kj  value  is  then  compared  with  Kj^  to  determine  the  load  (or  stress)  at 
initiation  of  crack  extension.  Such  a calculation  is  repeated  for  all  possible 
combinations  of  flaw  size.  The  result  is  a series  of  curves  for  a/a^  versus 
temperature  with  crack  depth  as  a parameter.  Figure  4 (21)  illustrates  two  sets 
of  initiation  curves  for  quasistatic  and  dynamic  loading  cases.  A part- through- 
wall  crack  is  assumed  in  Figure  4.  A separate  set  of  curves  should  be  developed 
for  another  flaw  geometry,  e.g.,  throughwall  crack.  It  should  be  noted  that  a 
set  of  such  curves  can  be  developed  only  for  the  specified  Kj^  at  a specific 
temperature . 

Discussion : 

The  fracture  mechanics  solution  used  in  the  development  of  crack  initiation 
diagram  is  based  on  a flaw  in  a flat  plate  (20,  21).  Such  a solution  ignores 
consideration  of  tank  car  cylindrical  geometry.  Further,  it  does  not  make  a 
distinction  whether  the  flav/  is  oriented  in  the  circumferential  or  axial  direc- 
tion. For  an  accurate  analysis,  solution  for  part- throughwall  flaw  in  a cylin- 
drical geometry  should  be  used.  Further,  separate  solutions  for  axial  and 
circumferential  flaws  should  be  used.  Finally,  where  appropriate,  an  elastic - 
plastic  fracture  parameter,  e.g.,  J-integral  parameter,  should  be  used. 

The  YC  point  on  the  crack  arrest  diagram  is  clearly  well  into  the  elastic -plastic 
fracture  region.  Consequently,  an  elastic -plastic  fracture  parameter  should  be 
used  to  define  this  point.  Further,  the  L-YC  curve  is  assumed  as  a straight 
line.  Accurate  L-YC  curves  should  be  developed  using  EPFM  method. 

A large  number  of  studies  have  shown  that  the  use  of  Kj  solution  and  Kj^  instead 
of  Kj.  can  lead  to  very  conservative  analysis,  particularly  when  ductile  fracture 
is  expected.  Reference  22,  among  others,  discussed  the  advantages  of  an  EPFM 
analysis.  It  is  well  known  that  fracture  may  not  occur  at  the  initiation  of 
crack  extension  when  toughness  increases  appreciably  with  crack  extension. 
Therefore,  improvements  can  be  made  to  initiation  curves  by  using  resistance 
curves  in  fracture  analysis.  A resistance  curve  (Kj-Aa  or  J-Aa)  based  analysis 
can  take  credit  for  additional  material  resistance  to  crack  extension. 

The  J-  integral  tearing  modulus  (J-T)  method,  which  applies  to  both  the  LEFM  and 
EPFM  fracture  regimes,  has  been  used  to  determine  the  load  at  which  unstable 
crack  extension  would  occur.  The  J-T  method  is  well  suited  to  examine  the 
possibility  of  fracture  at  any  point  at  and  beyond  initiation  of  crack  extension. 

The  analyses  performed  in  the  Evaluation  of  Railroad  Tank  Cars  section  illustrate 
how  an  EPFM  method  can  be  used  to  develop  critical  crack  size  curves.  These 
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curves  would  replace  the  initiation  curves  of  Figure  4.  The  section  XI  of  the 
ASME  Code  (23)  contains  analysis  methods,  procedures,  and  easy-to-use  flaw  eval- 
uation tables  for  carbon  steel  and  stainless  steel  piping  materials  and  associ- 
ated weldments  Refs.  24-29.  The  flaw  evaluation  procedures  for  carbon  steel 
piping  Refs.  27,  28,  29  include  screening  criteria  that  define  whether  a LEFM  or 
an  EPFM  analysis  is  appropriate  before  using  flaw  evaluation  tables  or  performing 
application  specific  analysis.  The  basic  elements  of  these  procedures  are  now 
being  considered  by  the  petrochemical  industry.  The  basic  elements  of  ASME 
procedures  can  be  developed  further  for  application  to  railroad  tank  cars. 

EVALUATION  OF  RAILROAD  TANK  CARS 


The  analyses  performed  in  this  work  are  based  on  the  J- integral  elastic -plastic 
fracture  parameter.  The  J based  analyses  are  superior  to  those  based  on  linear 
elastic  fracture  mechanics  concept  because  it  accounts  for  crack  tip  plasticity 
consistent  with  the  material  stress  - strain  curve  and  allows  assessment  of  the 
stability  of  crack  extension.  The  fracture  mechanics  methods  and  analyses  pre- 
viously used  for  the  tank  car  evaluation  were  based  on  LEFM.  Further,  the  solu- 
tions used  for  the  tank  car  were  for  the  simplified  case  of  a flaw  in  a flat 
plate.  Consequently,  the  effects  of  tank  shell  geometry  were  ignored  in  the 
analyses.  The  analyses  performed  in  this  section  use  appropriate  J- integral 
elastic -plastic  solutions  for  the  tank  car  geometry,  thus  allowing  accurate 
evaluation. 

Circumferential  Throughwall  Crack  Under  Axial  Tension : 

Stress  Intensity  Factor  Solution; 

The  mode  I stress  intensity  factor  for  the  tank  car  shell  geometry  with  radius -to 
thickness  ratio  (R/t)  of  100  is  given  below  (30): 

K,  = (a^)-(7rR^)'^^-F^(0)  (5) 

where:  >=  P/(27rRt)  and 

F^(^)  =1-1-  4.93- [5.3303- (^/7r)l-5  -i-  18 .773  ■ (0/n)^-^^] 

P is  the  applied  axial  load;  R,  t,  and  & are  the  shell  mean  radius,  wall 
thickness,  and  crack  half -angle,  respectively.  Figure  5 shows  the  through  wall 
crack  geometry.  This  solution  applies  to  0 < S/n  < 0.15  and  R/t  = 100.  The 

solution  has  the  following  accuracy:  (i)  better  than  2%  for  & /n  < 0.05  as  com- 

pared to  the  results  in  Ref.  9,  and  (ii),  the  solution  is  increasingly  con- 
servative for  0 /n  > 0.05. 

J- Integral  Solution: 

The  J-integral  solution  for  R/t  = 100  is  given  below  (30). 

J = (l/E)a^2(^R^^)p^2(^^)  + . ^) (B^(^))'"- (P/P^(5))'^1  (6) 
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where : 


Sq-  e [1  + (0.5F^(^)-)-{(n  - l)/(n  + l)]-{a^/a^)^  /(I  + (P/P^(^))2}] 

Pp(^)  = 2a^Rt[7r  -6-2  arcsin(0.5  sin^)] 

B^(^)  = [1  + 4.93  {5.3303(^/7r)l-5  + 18 . 773  ( 0/7r)^-2^  } ] / 

[1  + 1.1495  {5.3303(^/7r)l-5  + 18 . 773  ( ^/7r)^-2^  ) ] 

m = (n)0-5 

This  result  has  the  correct  elastic  limit,  shown  in  Eq.(5),  because  Oq  = 6 and 
Q = o,  when  n = 1.  H.,  depends  upon  6 /it,  n,  and  R/t.  In  the  above  solution,  H., 
is  to  be  used  for  the  R/t  '=  20.  The  solution  assumes  plane  stress  crack- tip 
condition;  the  plane  stress  condition  is  appropriate  for  the  tank  car  shell 
geometry.  The  R,  t,  and  6 are  the  tank  mean  radius,  wall  thickness,  and  crack 
half -angle,  respectively.  The  a,  a and  n are  constants  in  the  Ramberg- 
Osgood  stress-strain  relation: 

e/^o  ^ ^ OLia/o^)"^ . (7) 

The  and  are  the  reference  stress  and  reference  strain,  respectively,  and 
are  related  to  each  other  by  the  relation  o^/e^  = E.  The  is  0.2  percent 
offset  yield  strength  and  E is  the  elastic  modulus.  The  H.|  values  are  tabulated 
in  Ref.  5 for  R/t  = 20  and  r = 1,  2,  3,  5,  and  7.  Linear  interpolation  between 
specified  values  of  6 /n  and  n in  tables  is  permitted.  The  J solution  given  above 
is  applicable  to  n < 10. 

Loading  Cases: 

Several  loading  cases  were  considered  for  the  determination  of  critical  crack 
size.  These  included  loads  anticipated  or  implied  in  the  service,  test,  certain 
accident  condition,  and  design  condition  for  the  railroad  tank  car.  These 
loadings  are  defined  in  the  Manual  of  Standards  & Recommended  Practices  (31). 
The  loading  categories  considered  in  this  report  are:  (1)  impact  force,  (2)  test 
procedure,  (3)  start- to-discharge  pressure,  (4)  burst  pressure,  and  (5)  vapor 
tight  minimum  pressure.  Table  1 shows  the  loads  used  for  each  category.  These 
values  were  obtained  from  reference  (31).  Except  for  the  impact  force  case,  the 
load  was  defined  in  terms  of  the  internal  pressure.  The  axial  load  was  calcu- 
lated using  ttR^p  where  p is  the  internal  pressure  in  the  tank  car.  The  R and  t 
values  used  were  140  cm  (55  inch)  and  1.43  cm  (0.5625  inch),  respectively. 

Fracture  Mechanics  Analyses: 

The  fracture  mechanics  analyses  performed  in  this  report  are  based  on  the  J- 
integral  tearing  modulus  (J-T)  method.  This  approach  was  described  previously. 
The  method  compares  the  crack  driving  potential  (J)  to  the  material's  resistance 
(J-R)  curve.  If  the  calculated  J value  is  smaller  than  the  material  J value  at 
the  initiation  of  crack  extension  Jj  value,  then  crack  extension  is  implied. 
However,  this  crack  extension  can  occur  in  a stable  or  unstable  manner.  To 
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determine  this  possibility,  the  rate  increase  in  the  calculated  (applied)  J is 
compared  to  the  rate  increase  in  the  material's  resistance  to  crack  extension  at 
the  calculated  J value.  This  quantity  when  normalized  with  certain  combinations 
of  material  properties  is  called  the  tearing  modulus,  as  discussed  in  a previous 
section.  If  the  calculated  value  is  greater  than  that  for  the  material,  the 
crack  is  predicted  to  extend  in  an  unstable  manner  causing  fracture  of  the  tank 
car . 

The  J-T  method  is  applied  to  situations  where  ductile  tearing  is  anticipated. 
However,  when  ductile  crack  extension  is  not  anticipated  due  to  the  material 
behavior  at  the  temperature  of  evaluation,  the  method  is  still  applicable.  Here 
only  the  J values  are  compared  with  Jj  to  determine  the  critical  crack  size  for 
specified  load  or  critical  load  for  the  specified  crack  size. 

The  fracture  mechanics  analyses  were  performed  for  postulated  circumferential 
throughwall  cracks  in  the  railroad  tank  car  shell  region.  The  mean  radius  and 
the  wall  thickness  of  the  tank  car  used  were  140  cm  and  1.43  cm  (55.0  inch  and 
0.56  inch),  respectively.  The  analyses  were  performed  for  the  AAR  TC128  grade 
B normalized  material  at  service  temperatures  of  -40  °C  and  22  °C.  Only  the 
initiation  J values  were  available  at  these  temperatures.  Therefore,  the 
critical  crack  size  calculations  were  performed  for  applied  J values  equal  to  the 
material's  initiation  (Jj)  value. 

The  material  J values  at  -40  °C  and  22  °C  were  41  and  121  KJoules/M^  (240  and 
700  in-lb/in^),  respectively  (32).  The  J analyses  assumed  that  the  crack  was 
subjected  to  quasi-static  loading. 

The  J analyses  require  the  stress  - strain  curve  of  the  material  at  the  appropriate 
service  temperature.  In  particular,  this  requires  the  Ramberg-Osgood  (R-0) 
parameters  for  the  stress  - strain  curve  of  the  material.  The  values  of  R-0 
parameters  are  given  in  Table  2. 

The  procedure  for  determining  the  crack  size  at  the  initiation  of  crack  extension 
for  specified  axial  load  is  described  next.  The  applied  J- integral  is  calculated 
as  a function  of  crack  size  for  the  specified  load  and  tank  car  shell  geometry. 
The  applied  J is  calculated  using  the  solution  given  in  Section  entitled 
Circumferential  Throughwall  Crack  Under  Axial  Tension.  A large  number  of  crack 
lengths  were  postulated.  The  applied  J values  are  plotted  against  postulated 
crack  lengths  as  shown  in  Figure  6.  A line  representing  the  material  J.,  is  next 
drawn  on  such  a plot.  The  intersection  of  the  applied  J curves  with  the  J.,  line 
determines  the  critical  crack  size  for  the  specified  axial  load.  The  procedure 
is  the  same  for  any  other  type  of  load  or  flaw  geometry.  The  Ductile  Fracture 
Handbook,  Ref.  5,  gives  a sample  problem  and  the  step-by-step  solution  for  the 
bending  moment  loading  case . 

RESULTS 


Figure  6 shows  the  calculated  J values  for  five  axial  loading  cases  that  cover 
service,  test,  and  design  conditions.  These  results  are  for  -40  °C  service 
temperature.  As  mentioned  before,  the  calculations  were  performed  for  the  AAR 
TC128  grade  B normalized  steel.  The  critical  crack  size  is  determined  from 
Figure  6 by  obtaining  the  intersection  of  41  KJoules/M^  with  the  applied  J 
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curves.  The  critical  crack  lengths  are  summarized  in  Table  3.  The  calculated 
J values  for  service  temperature  of  22  °C  are  shown  in  Figures  7 and  8.  At  this 
temperature,  the  Jj  value  was  taken  as  121  KJoules/M^  (28).  This  is  a conserva- 
tive value  since  it  represents  J toughness  at  -18  °C.  The  critical  crack  size 
results  are  summarized  in  Table  3.  These  results  are  also  plotted  in  Figure  9. 

This  figure  can  be  used  to  determine  the  critical  crack  length  for  any  load  that 
is  bounded  by  the  loading  considered  herein. 

DISCUSSIONS.  CONCLUSIONS.  AND  RECOMMENDATIONS 

Fracture  mechanics  analyses  using  the  J- integral  as  fracture  parameter  were 
performed  for  two  service  temperatures  of  -40  °C  and  22  °C  and  five  axial  loading 
cases.  The  loading  cases  were:  (1)  coupler  impact  force,  (2)  test  pressure, 
(3)  start- to-discharge  pressure,  (4)  burst  pressure,  and  (5)  vapor  tight  minimum 
pressure.  Circumferential  throughwall  cracks  in  the  tank  car  shell  region  were 
postulated  to  determine  the  critical  crack  size. 

The  analyses  were  based  on  the  applied  J reaching  Jj,  the  material  toughness  at 
crack  extension.  Material  J- resistance  (J-R)  curves  could  have  been  used  in 
fracture  analyses,  but  these  were  not  available.  However,  the  J-R  curve  data  is 
relevant  above  NDT  temperature.  Consequently,  the  results  for  -40  °C  are  accu- 
rate. When  J-R  curve  is  used  in  the  analysis,  the  critical  crack  size  obtained 
is  larger  than  those  when  J is  used.  Consequently,  the  results  for  22  °C  are 
conservative.  The  J analyses  assumed  that  the  crack  was  subjected  to  quasi - 
static  loading. 

The  analyses  indicate  that  the  critical  crack  lengths  at  22  °C  are  larger  by 
factors  of  1.5  to  2 than  those  for  -40  °C.  The  critical  crack  lengths  for  dis- 
charge pressure  were  17.1  and  29.2  inch  at  service  temperatures  of  -40  °C  and 
22  °C.  These  crack  lengths  are  approximately  30  to  50  times  the  tank  car  shell 
wall  thickness.  The  critical  crack  lengths  for  the  coupler  impact  force  are 
approximately  1.7  times  larger  than  those  for  the  discharge  pressure.  The 
critical  crack  length  is  determined  to  be  the  largest  for  the  impact  force  case. 
The  results  also  show  that  even  at  burst  pressure  small  throughwall  cracks  (four 
times  tank  shell  thickness)  can  be  tolerated. 

In  a leak-before -break  scenario,  when  a part- throughwall  flaw  extends  through- 
the-wall  the  resulting  throughwall  crack  must  be  smaller  than  the  critical 
throughwall  crack  length.  While  analyses  for  part  - throughwall  flaw  were  not 
performed,  it  can  be  argued  that  a propagating  surface  flaw  having  length- to- 
depth  ratio  (aspect  ratio)  less  than  25  will  produce  leak-before -break  at  service 
temperatures  of  -40  °C. 

It  is  recommended  that  service  failure  data  should  be  examined  to  determine 
whether  part  throughwall  flaws  with  aspect  ratio  greater  than  25  can  exist  in  the 
tank  car  shell  region.  Fracture  mechanics  analyses  for  part  - throughwall  flaws 
are  recommended  to  determine  the  flaw  shape  and  size  that  would  produce  leak- 
before-break  condition  for  the  tank  car.  A safety  margin  should  be  applied  on 
the  predicted  crack  length  to  account  for  the  variation  in  J toughness  at  the 
evaluation  temperature . 
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TABLE  1.  Loads  Used  in  Fracture  Mechanics  Analyses. 


Loading  Category 

Pressure 

MPa  (psi) 

Calculated 

MN 

Axial  Load 
(lb-E6) 

Impact  force 

5.8 

1.3 

Test  pressure 

2.3 

340 

14.3 

3.2 

Start- to-discharge  pressure 

1.8 

255 

10.8 

2.4 

Burst  pressure 

5.9 

850 

36.2 

8.1 

Vapor  tight  minimum  pressure 

1.4 

204 

8.5 

1.9 

TABLE  2.  Material  Stress  - Strain  Parameters 


Temperature 


Parameter 

-40  “C 

22  °C 

a 

0.83 

0.91 

n 

7.93 

7.74 

MPa(ksi) 

488.5  (70.9) 

446.5  (64.8) 

MPa(ksi) 

608.4  (88.3) 

560.2  (81.3) 

E GPa(ksi) 

206.7  (30,000) 

206.7  (30,000) 

TABLE  3.  Critical  Crack  Size  for  AAR  TC128  Grade  B Tank 
Car  Shell,  Normalized  Steel. 


Loading  Category 

Critical 
-40  ° 

Throughwall 

C 

Crack 

22 

Length 

°C 

cm 

inch 

cm 

inch 

Impact  force 

85.3 

33.6 

127.3 

50.1 

Test  pressure 

31.5 

12.4 

55.9 

22.0 

Start- to-discharge  pressure 

43.4 

17.1 

74.2 

29.2 

Burst  pressure 

6.4 

2.5 

14.7 

5.8 

Vapor  tight  minimum  pressure 

55.9 

22.0 

90.9 

35.8 
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Figure  1.  Illustrations  of  the  J-T  Procedure. 
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Figure  2.  Illustration  Showing  a Curve, 
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NOMINAL  STRESS  RELATIVE  TO  U, 
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Figure  3.  AAR  Design  Reference  Graphs  for  L to  YC  Temperature  Transition 
Curves . 


16 


STRESS  RELATIVE  TO  Y.S. 
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Figure  4.  Illustration  of  AAR  Fracture  Analysis  Graph. 
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Figure  5.  Illustration  of  A Circumferential  Throughwall  Crack  and  Angular 
Relationships . 
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Figure  6.  Applied  J,  at  -40  C,  Versus  Crack  Length  for  Several  Axial  Loading 
Cases . 
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Figure  7.  Applied  J,  at  22  C,  Versus  Crack  Length  for  Several  Axial  Loading 
Cases . 
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J,  kJ/m''2 


Figure  8.  Applied  J,  at  22  C,  Versus  Crack  Length  for  P-  35.932MN 
(8,077,820  lb). 
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Critical  Crack  Length,  mm 


Axial  Load,  MN 


Figure  9,  Critical  Throughwall  Crack  Length  versus  load  for  normalized  AAR 
TC128  grade  B Tank  Car  Shell  Steel.  The  Curves  are  for  a 
Circumferential  Crack  in  the  Shell  Plate. 
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